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^^^V^^T^l^^^ -hing (or a wedge/ 
geometries and decomposes such into logifal^.t ^ *? app,ica b |e 
boundmg the decomposed pieces. It tien erTrce S Trie co k° Umes (b ^> ^im barrel caps 
mtersecting. and matching aliboundmg ball can^o^H C ° mpat,b ""y cons ™ by sprinting, 
throughout the geometry. 8 C3PS '° rap,d,y & nerate a well formed cont.nuous mesh 

1.0 Introduction 

» * » * m- W - {m f or aut0 . 
proposed solutions contains both positive and neS Zit^J" ^ d,fficult P^lems. each of the 
However powerful the tool's claimVare individual coun er , ? " ^ techni 1 ue ^P^d. 
definitely superior can be easily generic ^ HoSt Tr ? mP ' CS ° f ge ° metrieS where another tool is 
assortment may be able to realistically ZidTcZ'oflZfTT * int ° a ^.box 

or a tool manager) can determine a-priori the * the tools tn -^lves 

descnbed below, offers a unique automation JSS^^JST' " ^ ^ t0 ° l 

-^r^ today, society has a 

barrel is logicaNy a single axis CSSfc " S 

rows of elements) and on the ends by "caps" and " Su b C aos ! fto^u n $ ° f regUlar ^lateral 

pa-red faces). Although barrels are always logically s.nl T * eqU ' Va,em paired faces ° r ** of 

Thts technique enforces mesh compatibly bf en IIZLT f° metrica "y ™ arbitrary, 

barre cap using projection ( P 2S^^ST rT"" be,Wee " *" SUbC3pS 0n «* 
matchmg. It uses the barrel ribs to insure a validTan^^ t0 aCC ° mplish an opti ™' 

along the projection direction. Element connecSuXZ^J $ f ?\ em ' y generated interior nod « 
and thus actual element generation is extremely rapid "'^ determined b * the ca P meshes 

AswH.be described in this paper, the cooper too. has the following positive attributes. 
Complete automation for applicable geometries. 

• High element quality. 

• Boundary sensitivity (best elements along the boundary). 

• Orientation insensitivity (same mesh for transformed geometries). 
° Speed. 

i^W*:^^ - This includes a wide 

at,on of wedge elements (5 faced triangular pri m) 21 r'2 h " eqUa " y 3pp,icab,e t0 the S e "«- 

>ng triangle and/or mixed triang.e/quad^.ateral cap "th™ W * ge/heXahedral mesh « by si mp .y allow- 




2.0 Cooper Process 

The process of meshing using the cooper tool couples geometric evaluation, virtual decomposition, and pro- 

SSTri! rC hT S C r "f 011 ° f a gCneraI pUrp ° Se ' r0bust face meshin S tOGl ^ ^ the paving algo- 
rithm [1], and the use of a predictable regular mesh generation technique called submapping [9] and [ 10] 
This process involves the following steps: J ' 

1 . Classification of applicable volumes and sorting of faces into sides and caps. 

2. Meshing of all side faces using submapping. 

3. Formation and separation of barrel volumes. 

4. Iteratively matching and imprinting all subcaps to produce stacks of caps. 

5. Meshing of cap stacks. 

6. Projection of the barrel volume interior nodes. 

7. Formation of volume elements in each cap stack. 



3.0 Volume Classification 

In general terms, a vohime can be automatically meshed using the cooper tool if a logical decomposition can 
be found which would admit coupled single axis projections. For example, the volume in Figure 2 could be 
logically divided into 2 projections, each of which is a single axis projection. The cooper tool currently 
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requires ail shared "virtual" face* (th* f ^ 

between such projections to be ,„ ^T^nt vo.u^T " ? p""*' *°»^*«> occurred) 
» correctly classified by the cooper tool whi^'voh^ST T" F ' gUre 2a mee,s this c "«™ and 
detenmned by the cooper tool. To mesh the volume slown f Ficure ^ 33 ^ n °' " a * ai " «««'y 
ally decompose the original volume into piecesTacn of which W ° U ' d be ret > uired '° ™"u 

decomposuion is shown in Figure 3b and' Fig re 3c temnZ^™** ^ "* COO «* r «»•■ Such a 
natural next step in this work. g The ,,ft,n * of th,s com ™n projection restriction is 






a- Original Volum e h v - 

Hgure 2. Vbh^ ^k.u,. ..■„„ „ decomposition I Virtual d r-j, - 

face (dark shaded). ""*« t00 '- ^ P^ecuon direction is cons.stent along the slid vJ^T 






a - Qriginai volum p 1 . A 
Hgw. J. Ucron.posi.icn „/ a volume aSggg ' ' Acal .too,^ . , 

•nd 1 hecon, t . i „a, l on of.,, adjac* mash bo Jb^S i 

""""appable face, which ca, a. teSkSJ^I L T To ,llu!, ™« **• « 




not submappable. Determ ining a submappable face is accomplished bya simple calculation based on topol- 



a. Face geometry 



=^ =S faces, each of wh.ch ,s ^^T^l^^ is also 




_ I 1 

a. Face geometry 



ogy and angle designations and does not require a meshing attempt. 
A cap face is simply defined as any face which is not a side face. 




a. Cap faces front view 




SdX C,aSSmCa,1 ° n ° f 0,6 faCCS ° f *" exam «" e v ° lu ™ '"to »* faces^fgCshS Hp faces (dark 

iTme^ned C 3 r v e ;, a " SidC ^ - ^ 31 the requested densit y usi "g submapping. If cap faces 
are unmeshed, they mn.ally remam so to aUow greater flexibility during the batfel foJL'.nd meshing 



The Cooper Tool 



August 27. 1996 



4 



4.0 Barrel Formation 



ir s z^TL^:^ dihedrai •» * — ~ 

caps, an upper and a lower, and a set of exterior ribs ^ A ^ * defi " ed by two barrel 

side of the barrel. Due to tecoJS^^«£^ ° f * Which wi " P r °g™ along the out- 
logically down a barters side. Since eac InS 2^e 2* '* ^ theSe ribS Ca " then 
nodes, this process involves simply choosin S v „^" "connected logically to only four other 
p-^w^ 





— — a - Single rib extension 

n^m^ rCaCheS M " * b Nation ceases. A barrel is guaranteed to 

terminated into one or more £5 aps wh ich a^Tot h T C " P f "* <■* barrel has 

be reformed to contain subcaps S^^S^^^T^^^^^^^*^ 



a. Original volume 



b. Side face mesh 



c. First barrel 






p d a C °-"bi°edsubcap e. Second b arrel 

Figure 8. Formation of the barrels in a cooperable volume. 



f. Third barrel 
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5.0 Imprinting Subcaps 



As described above, each barrel will have an upper and a lower ran m. • 

» "hich logically connect exterior K£l nZl rt JeS ' gnat '° n ° f which ca P is flower 

rior barrel ribs. However, for all nontrivial Z^Z^T^Io Z^T " maintained With the «»* 
as shown in Figure 9. These subcaps rarely match an even when 1 , ^ C ° mp ° Sed ° f Subca P s 

matchmg operator A is denned for each barrel B such that a - B \ < a , 

sucn tnat a - A b (read subcap a matches 6 on 

barrel B ) and inversely b = *A 




^^^^^^^^^^^^^^ 

progresses these rib connections will include interior boundaries as well ow ever, as matching 

Th, to be ta prt nKd * fol classifie<i agains , ,„ ^ on ^ ^ ^ ^ bmi Lar , 

ITTT ' ' ^ " l0 ° P ' ° f ^ <"> * ' - — * Flgu* ,0, al 
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connected through barrel ribs as shown in Figure 10c. A subcap boundary loop can be 



divided into two sets: 




a. Schematic Upper Ca 




b. Schematic Lower Ca 




V which contains the portions of the loop which are attached to barre, ribs and V which contains the per- 
t,0nS ° ftheSUbCa P ,0 ° PS Whkh « " 0t ^ <° «*«. ribs as shown in Figure ... T ab can then be 





— ,, a - BMTelc "p b. Connected bounda ry 

Figure 1 1 Division of a barrel cap loop i " ~~ ~~ "~~ ~ - 




Unconnected Boundary 

barre.. The cap boundary shown ^tSS^SZ "odes to exist.ng nbs of the " 



> (b) and nodes without (c). 

defined as the subset of Y fl that is connected to T b through the existing ribs of the barre. as shown in 
Figure .2, Inversely, T ba is the portion of \ which is connected to T fl through the same ribs as shown 
■n F.gure .2b. T ab may be continuous or discontinuous, and may span multiple .oops of one/both of the 
caps. The characteristics of T ab determine the type of matching and/or imprinting that wil. occur 
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Upon _ successful completion of subcap imprinting, each subcap of the barrel will be paired with an eauiva 
lent subcap on the opposite end of the barrel, i.e. q 



(a= A b )V(a) 
(b = \)V(b) 

l h lrz e ^t zt*: p h oin ; e : on each cap/subcap - A cap/subcap may ^ « «* 

Ddrrei caps oecause it may be attached at most to two barrels. 

5.1 Subcap Matching 

If T ab is continuous on each loop i of cap*, {Y ab = rJ}Vi and {T ba = r' } V/ . then the two caps 
are one-to-one matches, and the barrel match pointers are set between the two caps a and b . If there are no 
subcaps on euher end of the barrel, then the caps by definition meet this criteria. 

5.2 Subcap Subtraction 

A subcap subtraction must occur when the projection of a subcap is completely within an opposing cap or 
subcap. More precisely, if T ab is continuous on one loop « of subcap a. T^tJ and 
{ r ac= 0} V(c * b) , then a subcap subtraction must be performed. This occurs when a subcap's bound- 
ary is connected along one portion of its boundary to only one other opposing cap/subcap A subcao is suh 
tracted by projectmg the "non-shared" boundary portion onto the opposmg ttn siting t 

opposing cap/subcap into two appropriate caps. Let t ab be the nonshared boundary portion of rj . or 



r ab = r'-r 



ab 



Y ab is then projected onto cap b to form a new boundary segment t ba . Two new caps, b and b' are now 
formed, where 

b={t ba + r ba ) 

The projected boundary segment, f ba is then connected to the original non-shared boundary portion f b 
by .ntenor nbs. These ribs contain only the two end nodes, and not any intermediate nodes along their length 
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through the barre,. Figure .3 shows one such subtraction. The new, y formed cap h is now , ma(ch w 
projected cap a and the matching barrel pointers are se, accordingly. 



o = V 






a. Cap fl to be imprin ted b. Opposing cap £ 

This example is for a rather simp.istic case, but the same principles hold for more complex cases such as 
when T ab spans multiple boundaries I of cap b or when / > 1 and ,h. «„. i • 

-P t/ or wnen / > I and the overlapping portion of the bound- 
ary is a complete loop of cap b or 



T ba = r/;0 e (2,/)) 



53 Subcap Intersections 



A subcap is said to intersect an opposing subcap if the projection is on.y partially in the opposing cap More 
preasely^foneloop^f^^ 

through ribs to two subcaps b and c such that the inverse of P a(b/c W is defined as 

r (b/c)a= ( r ba + r ca) 



then a subcap intersection exists and must be resolved. This case is shown in Figure 14. Let t a(b/c) be the 






a. Cap being imprinted b. Opposing interse ction caps c. Shared bounda ry portion 

F.gure 14. lntersect.on subcap schematic showing the intersection of barrel subcap a into opposing subcaps b and 



nonshared boundary portion of r ' , or 

From the definition of barrels, it can be shown that there will always exists an equivalent intersection case 
v.ewtng the barrel from the opposing end using either cap b or cap c . In other words, assuming we use sub- 
cap b , a boundary segment T b{a/d) will exist where subcap d is on the same barrel end cap as subcap a 
as shown in Figure 15. 









|T Cap d J 


ICap £>\ J 
















a. Opposing intersection caps 

Figure IS Inf^rcArfinn ciiiv>ork n n U«. 


b. Cap being imprinted 


c. Shared boundary portion | 



d. 



In either case, there is one common cap, cap a in this case, which forms the basis of an intersection. 
?a(b/c) is then projected onto caps b/d and the closest intersection node p is found. This intersection 
node must meet the even parity constraint if generating an all-quadrilateral/hexahedral mesh. Let Y b be 
the boundary segment connecting Y ba to the intersection point p , and let t ap be the projected portion of 

r a (b/c) also connecting r ba to the intersection point p as shown in Figure 16a. A new cap a is then 
defined as 



* = r ta + r bp + t ap 
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£EE523 r *" " ^ reVCrSe direCti0 " 35 Sh °- in , 6 and a second new cap 



ap 



n« pc„ io „ r„, »<,r v m to respeclively connccled by inlei . or , bs m r >nd ^ 
*» ln Figure lfc „„ lhe malches towra , fc5 , wo cips m m ^ _ ^ . ^ tf ^ 



a= A a . 






ately subtracted from caps a and & resnertivpiv Th»c a « u 
5.4 Subcap Seaming 



that nodes on edges of the externa, boundary are not combined by checking rib 



connections throughout the 




c - Imprint with a crack 

Figure 17. Seaming of projected subcaps 




d. Seaming of the crack 



5.5 Self-Intersecting Subcaps 

the two boundary portions that were combined to form cap ft-.r; and (T ' r 1 
intersect each other as shown in Fieure 18 b ba ' may overla P or 

a singie cap. the cap is said ^SZ^^S^ST *™ * Ending ,oop of 

of intersections, n , must be even The *™. ^ Cap " continu °^- the number 

the overling portions J^L^T^^ " ^ " + ' ~ ^ * 

dral mesh is desired to insure that the ^^^^2^.°" «■*» if » a"-hexahe- 

mation of the n + 1 new caps (cans A' h 77 ° f " 0deS - After the f °'" 

caps (caps b { b 2 and b 3 in Figure 18d) any completeiy interior caps ) 
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must be subtracted from adjacent caps (caps b and cap c in Figure ,8c) Th „■ 

are also shown in Figure I8d. } ' 11,656 adjusted ca P s *l and 




a- Upper subcap a being subtracted 





b-Lowersubcap* being imprinted 



on 




6.0 Meshing Barrels 
6.1 Forming Barrel Cap Stacks 

!g SStJ^S* - T -P <« - - mesh 
This mesh ,s then projected to each of the cansTnT Zl I"™" 3 ' USed for this m «hi"g Process 
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6.2 Generating Interior Barrel Nodes 

Since each ban-el now has a complete exterior mesh, this can be used to generate all the interior nodes of the 
mesh. Th,s ,s done us.ng a weighting of the same .east square residual method mentioned Son * Tut 
X,j be the least square residual projection of node i from the upper cap to layer ; . Similarly, let *,/ be 
the projection of node / from the lower cap to layer j . Then the location of node / on layer j is defined as 



f,3SS maintained in ° rder °" thC interiOT r ' bs as ^ « grated from the upper to the lower cap 
6.3 Volume Element Generation 

With all interior nodes generated, the formation of volume elements is a trivial ta* h«t ~rfi«« a ■ , 
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7.0 Examples and Evaluation 

^^ndfigure 19 through Figure 21 show 



examples of cooper tooi generated meshes. These 



examples 




a. From view 
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Several meshing densities for each of the examples are reported in Table I, but only the coarsest mesh is 
shown in the figures. The run times of the problems shown in this paper at the varying densities have been 
calculated on an HP9000 machine and are shown in Table I. The method is linear with the number of ele- 
ments, and roughly constant at over 100 elements per second. The slower times were mainly attributable to 
longer surface meshing times for nurb side faces. 



TABLE 1 . Cooper Tool Performance 



Figure 


Density 


Elements 


Elements/Second 


1/6 


0.35 


2,623 


133.7 


1/6 


0.20 


14,123 


126.5 


8 


1.00 


1,322 


142.9 


8 


0.70 


6,328 


138.3 


19 


0.70 


948 


81.0 


19 


0.30 


12,244 


101.84 


20 


1.00 


789 


48.3 


20 


0.70 


2,383 


48.9 


21 


1.00 


2,137 


121.0 


21 


0.60 


8,652 


137.4 



8.0 Conclusions 

A new cooper tool has been presented which automates hexahedral meshing of a large class of geometries. 
The method decomposes a volume into a unique set of logical barrels which are then meshed using a single 
axis projection. The method is fast and robust. Element quality is assured through the generation of consis- 
tent surface meshes for projection. The cooper tool meshes are boundary sensitive (best mesh at the bound- 
ary) and orientation insensitive (same mesh for transformed geometries). 
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